The pre-inflationary vacuum in the cosmic microwave background 
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We consider the effects on the primordial power spectrum of a period of radiation-dominated ex- 
pansion prior to the inflationary era. If inflation lasts a total of only 60 e-folds or so, the boundary 
condition for quantum modes cannot be taken in the short-wavelength limit as in the standard 
perturbation calculation. Instead, the boundary condition is set by the vacuum state of the prior 
radiation-dominated epoch, which only corresponds to the inflationary vacuum state in the ultra- 
violet limit. This altered vacuum state results in a modulation of the inflationary power spectrum. 
We calculate the modification to a best-fit model from the WMAP3 data set, and find that power 
is suppressed at large scales. The modified power spectrum is favored only very weakly by the 
WMAP3 temperature and polarization data. 



INTRODUCTION 

The high degree of spatial flatness, isotropy and ho- 
mogeneity that we observe in the universe today is well 
explained by inflationary cosmology [![. In order to ob- 
tain the universe that we observe, the minimum number 
of e- foldings of inflationary expansion must lie somewhere 
between N = 46 and N = 60, but inflation could have 
lasted for significantly longer. The nature of the universe 
prior to inflation is unknown, largely due to the fact that 
inflation is very good at washing out initial conditions. 
The existence of attractor solutions to the inflationary 
equations of motion make obtaining information about 
the pre-inflationary universe virtually impossible if infla- 
tion continues for long enough. Implementations of infla- 
tion in string theory (e.g., Ref. 0) generically suffer from 
the "rj problem" , which means they must be fine-tuned to 
achieve sufficient inflation. It is therefore well-motivated 
from the point of view of string theory to consider models 
with a minimal number of e-folds. If inflation lasts just 
long enough to solve the horizon and flatness problems, 
it is quite possible that pre-inflationary cosmology will 
leave imprints on the present day universe. In particu- 
lar, pre-inflationary physics can lead to modulations in 
the primordial power spectrum, which may be observed 
in the large-scale temperature anisotropics of the cosmic 
microwave background (CMB). 

During inflation, cosmological perturbations had their 
start as tiny vacuum fluctuations which were then red- 
shifted to super-horizon scales by the rapidly expanding 
spacetime. Perturbations on the scale of the present-day 
universe left the inflationary horizon at around N = 60, 
having originated as small-scale vacuum fluctuations at 
some earlier time. In order to fully reconstruct the pri- 
mordial spectrum generated by inflation, we must be 
able to understand the evolution of quantum modes from 
their birth in the inflationary vacuum out to the extreme 
large-scale limit. If inflation lasts for many more than 60 
efolds, which is commonly assumed, then the inflation- 
ary model fully specifies the evolution on all observable 
scales. However, if inflation lasts for just the minimum 



number of e-foldings, then the initial conditions of the 
quantum fluctuations must be set in the pre-inflationary 
era. Information about this earlier time is therefore re- 
quired to uniquely calculate the inflationary power spec- 
trum. 

In this paper, we study the effects of a prior period 
of radiation-dominated expansion on the subsequent in- 
flationary epoch in models which yield the minimum 
amount of inflation, which we take to be N = 60 for 
definitcness. In such a circumstance, the initial quan- 
tum modes are born not in the inflationary vacuum, but 
in the vacuum appropriate to radiation-dominated ex- 
pansion. This has a profound effect on the largest-scale 
perturbations. The standard calculation of perturbations 
in inflation assumes a well-defined short- wavelength limit 
for quantum modes. However, a fluctuation that leaves 
the horizon at the very start of inflation would need to 
have 'started out' large, since no non-inflationary mech- 
anism could have brought it to such scales. Therefore, 
when we set initial conditions for the fluctuations at the 
start of inflation, we are forced to consider quantum fluc- 
tuations far from the short-wavelength limit. As we show, 
fluctuations that are born at large wavelength during the 
radiation-dominated era affect the power spectrum very 
differently than fluctuations born at small wavelength 
during the inflationary era, resulting in a strong suppres- 
sion of power on large scales. 

The transition from a radiation-dominated phase to a 
de Sitter universe was first studied by Vilenkin and Ford 
Starobinsky Q, and by Linde Q, and the produc- 
tion of gravity waves during such a transition has been 
studied by Sahni 0. These earlier works concentrated 
on the impact of thermal effects on the stability of the de 
Sitter phase transition. In contrast, we consider the case 
of a decoupled scalar field and assume an instantaneous 
phase transition. Not only does this simplify the analy- 
sis, but it elucidates the profound effect that the choice 
of vacuum alone has on the primordial power spectrum. 
We emphasize the role of the choice of vacuum in the 
modulation of the power spectrum, as opposed to the 
pre-inflationary dynamics of the inflaton field itself, as 
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considered by previous authors @, H, @] ■ 

This paper is organized as follows: in Section I we re- 
view the standard treatment of inflationary fluctuations 
in de Sitter space. In Section II we study the evolution 
of both inflaton fluctuations and curvature perturbations 
in the pre-inflationary radiation-dominated universe. We 
find that if the inflaton is slowly rolling then the two ef- 
fectively decouple, allowing us to treat the inflaton as a 
massless scalar during radiation domination. This allows 
us to construct an adiabatic vacuum that is exact at all 
wavelengths. In Section III, we calculate the modification 
to the best-fit WMAP3 power spectrum when we initial- 
ize the quantum modes in this vacuum. We calculate the 
likelihood of the modified spectra for varying amounts of 
inflation and find a spread in likelihoods | — 2Aln£| < 2, 
the best-fitting model yielding — 2Aln£ = 1. Despite the 
strong modification to the underlying power spectrum, 
the Cj T - and C/^-spectra are altered only modestly, 
however, the Cj^-spectra are more sensitive to this ef- 
fect. Therefore, power spectra with suppressed ampli- 
tude at large scales are favored only weakly by current 
data, but improved measurements of the polarization sig- 
nal may be able to detect this effect amidst the statistical 
scatter. In Section IV we present conclusions. 



VACUUM SELECTION IN DE SITTER SPACE 

In this section we briefly review the standard calcula- 
tion of inflationary fluctuations, in which the ultraviolet 
limit is assumed to exist for all quantum modes. 

During inflation, quantum fluctuations in the field or 
fields driving the expansion are quickly redshifted to 
scales far greater than the causal horizon, where they 
manifest themselves as classical curvature perturbations. 
Specializing to the case of a single scalar field (f> = <fio+5if, 
the background field, 0o, couples at linear order to the 
metric perturbation. It is therefore useful to introduce 
the gauge invariant Mukhanov potential, 



(1) 



where a is the scale factor of the universe, H is the Hubble 
parameter and primes denote derivatives with respect to 
conformal time. The field fluctuation, S<p, and potential, 
'J, are gauge invariant and defined as in Ref. 



84) = 5ip + (f)' (B-E'), 
* = V- (^) (B-E% 



(2) 
(3) 



where the perturbed metric is, 

ds 2 = a 2 (r) [(1 + 2(j))dT 2 

- 2B ii dx i dr - ((1 - 2?A)% + 2E tij )dx i dx : >] . (4) 



In comoving gauge, 5<f> = 0, and the metric perturbation 
measures the spatial curvature perturbation on comoving 
hypersurfaces, ip = TZ. In this case we have 



K = 



(5) 



where z = 4>'/H. The power spectrum of the curvature 
perturbation may then be written as a function of co- 
moving wave number k, 



k 3 u k 2 



(0) 



where the u k are the Fourier modes of the gauge- invariant 
potential satisfying the equation 



k 2 - 



u k = 0. 



(7) 



A non-zero correlation function CR- 2 .) means that there 
has been particle production due to the inflationary ex- 
pansion. In an expanding spacetime, the initial vacuum 
state associated with positive frequency quantum modes 
will in general be associated with a mixture of both pos- 
itive and negative frequency modes at later times. An 
observer in a state initially devoid of quanta will regis- 
ter a thermal bath of particles as the universe expands. 
When solving Eq. ([7|) the initial conditions on the modes 
correspond to a specific choice of vacuum. As an exam- 
ple of vacuum selection, consider the simple case of a 
free scalar field fluctuation evolving in a de Sitter back- 
ground, H = const. It is convenient to work with the 
variable y — k/aH, which is the ratio of the Hubble ra- 
dius to the physical wavelength of the perturbation. In 
a de Sitter background the conformal time is 



I 

so that y — —kr and Eq. |(7J) becomes 



2 d 2 u k 
dy 2 



+ (y 2 - 2)u k = 0. 



This has the general solution 



'3/2 



3/2^ 



(8) 



(9) 



(10) 



where H^ 1 ' 2 ' are Hankel functions of the first and second 
kind, respectively. Since de Sitter space is eternal, the 
initial conditions can be set in the infinite past, y — > oo. 
In this limit Eq. (fTO)) becomes 



uk(y) 



1 



2k 



de ty + c 2 e 



(11) 



In order that the mode reduce to the vacuum state 
which annihilates positive frequency excitations, we must 
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choose C2 = 0, corresponding to the Bunch-Davies vac- 
uum. Then Eq. (|10|) becomes 



u k{y) = — 



1 



f2k 



1 



- ) 

V 



(12) 



This solution starts off as a plane wave in the infinite 
past and at late times as y — > 0, the second term in 
parentheses comes to dominate and particle production 
occurs. 

In general, not all exact solutions to Eq. ([7]) will 
asymptote to plane waves at initial times. In this case, 
the best one can do is construct an adiabatic vacuum to 
some finite order in a parameter which characterizes the 
slowness of the expansion. The adiabatic vacuum can 
be thought of as the vacuum that best approximates the 
Minkowski vacuum in an expanding spacetime. In order 
for such an approximation to be sensible, the background 
expansion must be slow relative to the frequency of the 
quantum fluctuation. From Eq. (J7|) with z" / z = C(r) 
and u> = \Jk 2 — C(r), this condition can be written 



<fln(C) 
dr 



<C LJ. 



(13) 



The vacuum obtained above clearly satisfies this con- 
dition (in the infinite past the cosmological expansion 
C(r) — > 0), and is equivalent to the lowest-order adia- 
batic vacuum. In de Sitter space, the adiabatic condition 
Eq. (fT3"|) is — kr ^> 0, which is not only satisfied in the 
infinite past, but at any later time for infinitely-large mo- 
mentum modes, k — > oo. Therefore, as long as we take 
the initial mode momenta sufficiently large, the adiabatic 
vacuum is a good approximation to the exact solution. 
This is to be expected, since at such small length scales 
the quantum modes do not feel the expansion of the uni- 
verse. However, modes that exit the horizon at the start 
of inflation will not have originated in the short wave- 
length limit and are not well described by the state Eq. 
(fTTjl . The 'initial state' of the quantum modes at the 
start of inflation must instead be determined by the pre- 
inflationary physics. In the next section, we consider the 
case of radiation-dominated (RD) expansion prior to the 
inflationary epoch. 



QUANTUM FLUCTUATIONS IN A RADIATION 
DOMINATED UNIVERSE 



In a radiation dominated universe, modes evolve under 
expansion from the long-wavelength limit y <C 1 to the 
short- wavelength limit t/ > 1. Quantum fluctuations on 
the scale of the inflationary horizon at the start of infla- 
tion will therefore not have originated in the Minkowski 
vacuum in the ultraviolet. We nevertheless expect quan- 
tum fluctuations to exist on all scales, including super- 
horizon, during the radiation-dominated phase 11| . Off- 
shell modes may be populated acausally as long as they 



do not couple to classical perturbations on superhorizon 
scales. These modes will be later "converted" to classical 
perturbations during inflation. In this section, we study 
the behavior of these fluctuations in the RD phase, their 
evolution determining the initial state for the ensuing in- 
flationary expansion. 

The dynamics of a decoupled scalar field that is mini- 
mally coupled to gravity is described by the action 



S 



-g d i x 



g^d^dA - V(</>) 



(14) 



Specializing to the case of a Robertson- Walker spacetime, 
variation with respect to the field <fi results in the Klein- 
Gordon equation, 



2aH(j)' + V^a 2 = 0, 



(15) 



where primes denote derivatives with respect to con- 
formal time. The equation of motion for the quan- 
tum fluctuations of the field are found by perturbing 
the field about the homogeneous background solution, 
(j) = 4>q + <50, and linearizing Eq. (|15p about this back- 
ground solution, 

8<t>" + 2aH6<j)' - V 2 ^ + ^a 2 (50-40'$' + 2^ a 2 $ = 0. 

(16) 

The background field couples at linear order to the gauge 
invariant metric perturbation, $, which is defined as in 
Ref. [ll|, 



$ = <f>+- [(B-E')a]' . 



(17) 



We assume the absence of anisotropic stress, so that 
^ = $. Because of this coupling, we must also study 
the evolution of the metric perturbation, Further- 
more, because the universe is radiation dominated at this 
time, there will generically be fluctuations in the radia- 
tion fluid as well. Perturbing Einstein's equations to first 
order, we obtain the equations of motion of the metric 
perturbation, 

V 2 $ - 3aff$' - 3a 2 H 2 $ = — (a 2 Sp r + 5 P4> ) , (18) 



in 



pi 



$" + 3ai?$' + 3a 2 H 2 $ + 2aiT$ 



4tt 



'pi 



(a 2 5p r + dpj,) , 
(19) 



where Sp r = 5p r /Z denote the pressure and density of 
the radiation fluid, and for the scalar field we have 



= -0' 2 $ + <t>'5<f>' + V> 2 



(20) 
(21) 



In this paper, we do not adopt a detailed model of pre- 
inflationary physics and so the exact nature of the curva- 
ture perturbation is unknown. However, one generically 
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expects that both thermal fluctuations in the radiation 
fluid as well as inflaton vacuum fluctuations will gener- 
ate a curvature perturbation across a range of scales. In 
what follows, we consider the case where the full range 
of comoving scales relevant for CMB physics are ini- 
tially superhorizon during radiation domination. With 
kg=2 / ki = 2ooo ~ 1CP 3 and if we suppose that ^=2000 en- 
ters the horizon just below the Planck scale, then kg = 2 
enters the horizon at H ~ 1CP 6 topi. Since we are consid- 
ering a minimum amount of inflation, we require that the 
quadrupole be on horizon scales at the start of inflation, 
giving flinf = 10 _6 TOpi. This assumption simplifies the 
analysis somewhat, but relaxing this condition does not 
significantly affect the qualitative results. 

We first study the initialization of the thermal fluctu- 
ations. The fractional thermal density perturbation in a 
comoving volume during RD is given by [l3T | , 



5p r 

Pr 



16 

35' 



(22) 



where S is the total entropy within the comoving region. 
The entropy density is given by 



2^ 2 t3 



(23) 



where g* s denotes the number of effective relativistic de- 
grees of freedom, and the temperature T can be written, 



T 



1 



1/4 



9 



Hmpi 
1.66 ' 



(24) 




Using these expressions with g* s sa the initial am- 
plitude of a thermal fluctuation that evolves with the 
expansion to a physical scale kj a is 



(25) 



where y = k/aH. The physical scale of interest is 
that attained by the fluctuation by the start of inflation, 
2/inf = k/a m iH ln {. Furthermore, we know that fluctua- 
tions on the scale of the quadrupole were on the scale of 
the horizon when inflation began, ki=2 = Qinf^inf- This 
gives 



te=2 



3/2 



Hj n f 

m P i 



3/4 



(26) 



It is important to note that although Eq. (j2"B"|) is to be 
evaluated at the beginning of inflation, it gives the initial 
amplitude of the thermal fluctuations on each comoving 
scale. Superhorizon modes should be initialized at hori- 
zon crossing, since it is unlikely that the universe will 
have been in thermal equilibrium on superhorizon scales 
prior to inflation. 



In order to determine the initial values of the inflaton 
fluctuations, one must make an assumption about the 
pre-inflationary vacuum state. If we suppose that the 
fluctuations originate in the Minkowski 'in-vacuum', a 
choice we will soon motivate, then they can be canon- 
ically normalized, S(f>i cx (oiVzfc) - 1 - As mentioned, 
although these fluctuations may exist on superhorizon 
scales, we do not expect them to couple to the curvature 
acausally. Like the thermal fluctuations, we therefore 
initialize them at horizon crossing. 

We are now ready to calculate the curvature perturba- 
tion generated by the combined effect of the thermal and 
vacuum fluctuations. Consider the i-0 perturbed Einstein 
equation, 

($' + aH<&) , = {Su + ^'S^i), (27) 

TOpi 

where Su denotes the velocity perturbation of the radi- 
ation fluid. Since thermal fluctuations are not vecto- 
rial, we can neglect the velocity perturbation at early 
times. Substituting Eq. (|2T)l into Eq. (fT5)l and moving 
to Fourier space yields the equation, 



47T 



47T 



- k' 



uipi V'^pi 

[a 2 5prk + 4 



54>' k + (V^a 2 + 3aHtf) 5<j) k ] .(28) 

The thermal and inflaton fluctuations enter this equation 
as sources proportional to k 2 5Tk and (V^a 2 +3aH(j)')6<fik, 
respectively. With <j>' cx H ln i^Je and cx H 2 nf ^/e, the 
relative initial amplitudes of the sources may then be 
written, 

k? <hk_ _ mv\H t {k)8 Tk 

(U a 2 + 3aH<t/)6<l> k 




where Hi (fc) is the value of the Hubble parameter when 
mode k enters the horizon. We see that the initial 
thermal fluctuations easily dominate the initial vacuum 
fluctuations. For example, with H m { — 1CP 6 topi and 
e = 0.001, the above ratio 
remaining terms involving 



nf 

i 2 x 10 4 (fc/fc £=2 ) 3 / 2 . The 
and (f)' 2 , being proportional 



t° H 2 nt \Je/k and H 2 n{ e/k 2 , respectively, are also sup- 
pressed. We therefore expect the inflaton fluctuations to 
decouple from the evolution of the pre-inflationary cur- 
vature perturbation. In this circumstance, the curvature 
perturbation evolves as it would in a single-component 
radiation universe, with the decaying solution [ll| . 



$(x,t) 



kr ( kr \ 



sin 



kr 

71 



kr . I kr 



v'k-x 



(30) 
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Therefore, if the scalar and radiation perturbations are 
decoupled at horizon entry, they will remain decoupled 
as the modes evolve to subhorizon scales during the 
radiation-dominated evolution. 

We are now in a position to determine the influence 
of $ on 8<fi. We can gain insight into this question by 
comparing the magnitudes of the source terms in Eq. 
(fT6|) . Since $ is decaying rapidly during RD, it suffices 
to compare the initial values of the source terms. Since 
the thermal fluctuations are dominant, Eq. (|27|) reduces 
to the Poisson equation, 



4tt 2 
—a park- 



in 



pi 



The ratio of source terms is then 



2V>?$fcj ( Hi 



k 2 c 



inf 



'JTk 



Hi{k)J Sfa 



(31) 



(32) 



where we have neglected the term involving V t ^, which is 
small during slow-roll. By substituting Sxk /54>k from Eq. 
(|2"9")l . we find that this ratio can be easily made <C 1. The 
term involving in Eq. (|I6p is proportional H ln iek$>i 
and can be dropped, with the result that the S(j)k are 
decoupled at linear order from the Newtonian potential 
& k created by the thermal perturbations in the radia- 
tion. The physical origin of this decoupling is easy to 
understand: during the period when the energy density 
of the universe is dominated by radiation, the Newtonian 
potential is likewise dominated by perturbations in the 
radiation field. However, in the limit that y^e ~ V' /V is 
small, the scalar is to a good approximation decoupled 
from the Newtonian potential at linear order and evolves 
as a free field in the radiation-dominated background. 
We show below that, despite the fact that the potential 
is dominated by thermal fluctuations prior to inflation, 
the curvature perturbation after inflation is dominated 
by fluctuations in the scalar field. This is again true as 
long as e is small. It is therefore self-consistent to neglect 
thermal radiation perturbations in the analysis which fol- 
lows, as long as the scalar field potential is sufficiently 
flat. 

With these simplifications, the equation of motion for 
the inflaton fluctuations Eq. (fTo]) becomes 



b'l + 2aH5(p' k 



< -r k 2 S(j) k + V^a 2 l 



■k = o. 



(33) 



Specializing to the case of radiation dominated expan- 
sion, a oc r, and redefining the field, u k — a8<j) kl this 
becomes 



u'l + (k 2 + a 2 V u ) u k = 0. 



(34) 



As mentioned, it is safe to neglect the term proportional 
to V i( p<p if the following inequality is satisfied: 



1 > 377, 



(35) 



where rj w m^V^ / (8wV) . This is valid during RD as 
long as r] does not vary significantly, a good approxima- 
tion during slow-roll. If we then change our time variable, 
y = kr, we obtain the expression, 



d 2 u k 
dy 2 



u k = 0. 



(36) 



This is the equation of motion of a massless scalar field in 
a radiation dominated background, and the main result 
of this section. This has the plane wave solution 



u k = cie 



-ly 



c 2 e 



(37) 



This demonstrates the well-known fact that there is no 
particle production during radiation-dominated expan- 
sion, since the quantum fluctuation never evolves out of 
its vacuum state. Since C(t) — 0, the adiabatic condi- 
tion Eq. (fT3|) is satisfied for all time and for all k. The 
vacuum is therefore well defined for all momenta and it 
remains to determine the coefficients c\ and ci . We know 
that at in the ultraviolet limit, y = kr — * oo, the modes 
should reduce to the Bunch-Davics form, 



u k (y) 



1 



2k 



(38) 



However, it is not clear what the initial state of the large 
scale fluctuations should be. From the field decomposi- 
tion, 



d 3 k / \ 
S0(x,r)= / ( 27r )3/2 { u k( T W + Mfc(-r)a T _ k j 



,ik-x 



(39) 

and the solution Eq. ([37]) . we obtain the alternative de- 
composition, 



5<f)(x, t) 



(40) 



Here we have introduced the Bogoliubov rotated ladder 
operators, 



6 k = cia k + Cjalk, 
= cjal k + c 2 a k . 



(41) 



A reasonable choice for an initial vacuum is 6 k |0) = 0, 
from which it follows that c 2 = for all k. This particular 
choice of vacuum corresponds to the state of minimum 
uncertainty, and has been used to study trans-Planckian 
effects [lj, [l5|, as well as the quantum decoherence of 
initial inhomogeneities [l6j]. Once we choose c 2 = 0, the 
value of C\ is determined from the commutation relation 
between the field fluctuation and its conjugate momen- 
tum, 

[<ty(x, t), tt(x', t')W = *£ 3 (x - x'), (42) 
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where it = a 2 S(f>' . With this choice of vacuum, we obtain 
the mode functions during RD, 

1 



Uk(y) 



'2k 



(43) 



Other choices for the initial state of the fluctuations are 
possible, and one ultimately must make an assumption 
about the pre-inflationary physics. This particular choice 
is pure positive frequency on all scales, and asymptotes 
to the Bunch-Davies limit in the ultraviolet. 

While we have determined that the curvature and in- 
flaton perturbations evolve independently during RD, we 
must not forget that the pre-inflationary curvature spec- 
trum may contribute to the post-inflationary spectrum 
generated by inflation. Since the pre-inflationary curva- 
ture perturbation decays (cf. Eq. (|30|) ). we expect a 
negligible contribution on small scales. However, on the 
largest scales, the pre-inflationary curvature perturbation 
can be approximated by 



P, 



k 3 

K,pre(fc) = ^l$fc*| 2 . 
Z7T 



(44) 



The inflation-produced curvature perturbation spectrum 
is P K (k) = H 2 /{izeml x ). With $ fel given by Eq. (31} 
and considering scales k ~ ki =2 , we obtain the ratio, 



■R-.prc 



(k) 



Pn{k) 



cx lOe 



k ~- k ,• 



mpi 



H 



inf 



1/2 



For flinf = 10~ 6 m P i, we find that e < 10~ 
trum produced by inflation is to dominate. 



(45) 
if the spec- 



THE POWER SPECTRUM AND ITS EFFECT ON 
THE CMB 

We are now in a position to determine the power spec- 
trum that results from such a short period of inflation. 
The idea is to set the initial conditions, Eq. (|43[) . for each 
mode at the start of inflation, and then evolve them with 
Eq. 0. As a concrete example, consider the exactly sol- 
uble case of power law inflation, for which a ~ t 1 ^ and 
e = const. In this background, z" / z = a" /a, and the 
general solution to Eq. (7]) is 

u k (kr) = V^fr? ( Cl ff«(-fcr) + c 2 ff( 2 )(-fcr)) , (46) 

with v — (3 — e)/(2 — e). The coefficients are determined 
by matching to the appropriate initial conditions, i.e. 
selecting a vacuum. The power spectrum is obtained by 
evaluating Eq. (6]) in the asymptotic limit, — fcr — > 0. In 
this limit, Eq. (j46|) becomes, 



Ufc(fcr) 



-kr- 



iV^krY v (-kT)(ci - ca), (47) 



where Y v is a Bessel function of the 2 nd kind. When c\ = 
(l/2)(yV/fc) and c 2 = 0, we recover the standard result 



obtained when the quantum modes are assumed to have 
originated in the Bunch-Davies vacuum. We therefore 
rewrite Eq. (|4"7} as 



u k (kr) 



-kr — >oc 



ubd(ci - c 2 ), 



(48) 



where ubd denotes the asymptotic form of the Bunch- 
Davies mode solution. The power spectrum follows: 



Pn(k) 



C2 



c x c 2 - CiC 2 ) 



(49) 



This expression gives the modified power spectrum in 
terms of the standard spectrum, Pj^ , and a vacuum 
dependent term proportional to the coefficients. This is 
valid for any vacuum choice, with c\ and c 2 generally 
being complicated functions of k. 

We now apply this approach to determine how the 
WMAP3 best-fit spectrum is modified if we suppose that 
inflation was preceded by a radiation-dominated epoch. 
We take the best-fit model with scala r sp ectral index 
n s = 0.951, no tensors and no running [l2|. We assume 
that the best-fit spectrum is generated by some inflation 
model with sufficient inflation to ensure that all observ- 
able scales originated in the Bunch-Davies limit of the 
inflationary vacuum. This is equivalent to specifying the 
observable n s at N = 60. As noted earlier, scales corre- 
sponding to the quadrupole can exit the horizon at any 
time between N = 46 and N — 60. Here we choose 
N = 60 but our results, Figures I-IV, are not sensitive 
to this specific choice. What is important is the amount 
of inflation occurring before the quadrupole leaves the 
horizon. We compare this result to the case of the same 
best-fit inflationary model, but with a pre-inflationary 
RD phase ending at N = 60 (Figure Q} . Instead of using 
the Bunch-Davies boundary condition for the modes, we 
use the RD vacuum solution, Eq. J43|) . as a boundary 
condition. The suppression of large-scale power is imme- 
diately evident in Figure [TJ The scales corresponding to 
these fluctuations are vacuum modes that were in exis- 
tence on scales of order the horizon size at the onset of 
inflation. As a result, their amplitudes are highly sup- 
pressed in relation to quantum modes that would have at- 
tained these scales as a result of inflation. Modes evolving 
during inflation undergo mode freezing as they cross out- 
side of the horizon, effectively locking-in their amplitudes 
at relatively large values. On sub-horizon scales, k > aH, 
the modified spectrum undergoes oscillations, rapidly ap- 
proaching the standard inflationary spectrum. While the 
inflationary Bunch-Davies vacuum only exists for modes 
k ^> aH, it has the same form as Eq. . 1 Therefore, 



1 There is an overall unobservable phase shift between the real 
and imaginary parts of owing to the fact that — kr — > kr in a 
radiation dominated phase. 
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FIG. 1: Best-fit WMAP3 spectrum calculated using standard 
inflationary initial conditions along with the same inflation 
model but with a radiation-dominated pre-inflationary phase. 



as k — > oo, the inflationary vacuum approaches the RD 
vacuum, and the mode solutions become identical. 

A similar spectrum was obtained in 0], in which the 
effects of an initially fast-rolling inflaton were studied. 
Burgess et al. [8( also obtain similar spectra arising from 
a hybrid model in which a rapidly oscillating auxiliary 
field leads to an era of matter domination prior to in- 
flation. We emphasize that we are encoding the effect of 
the pre-inflationary physics solely in the choice of vacuum 
state for the inflaton fluctuations. We do not address the 
separate questions of how the universe achieves the neces- 
sary level of homogeneity prior to the onset of inflation, 
or of how the universe makes the transition from the 
radiation-dominated to the inflationary phase (the for- 
mer could be achieved via a previous period of inflation, 
as in "open" inflation models.) In this sense, our analysis 
should be regarded as a physically motivated "toy" model 
which demonstrates the effect of the pre-inflationary vac- 
uum state on the primordial power spectrum. 

Using a modified version of CAMB [l7|, we are able to 
generate the CV-spectra for prior RD. In Figure 2, we 
plot the temperature spectrum. In each model consid- 
ered in Figure 2, scales corresponding to the quadrupole 
leave the horizon at N — 60. These models differ in 
the total amount of inflation they provide, as indicated 
in the figure. The lack of power at large-scales is ev- 
ident in the low-£ multipoles, but the spectra become 
virtually the same at around I ~ 11. The modification 
to the spectrum is strongly dependent on the duration 
of inflation, and the effect is completely washed out for 
models yielding more than 1 e-fold of inflation between 
the onset of inflation and the time when scales corre- 
sponding to the quadrupole leave the horizon (N > 61 
in the figure.) In Figures 3 and 4 we plot the polariza- 



FIG. 2: Cj -spectra for the best-fit WMAP3 model (solid 
black) together with models for which inflation is preceded 
by radiation dominated expansion. Models are plotted for 
varying durations of inflation. Note the suppression of power 
at large scales. 



tion auto-correlation (EE) and temperature-polarization 
cross-correlation (TE) spectra. Each model depicted in 
these figures lies within | — 2Aln£| < 2 of the best- 
fit model, and cannot be distinguished statistically with 
high confidence. The temperature spectra all lie within 
the cosmic variance envelope, making this effect difficult, 
if not impossible, to decisively resolve with future exper- 
iments. The EE spectra are the least affected by the 
pre-inflation phase. However, the TE cross-correlation 
spectrum is more sensitive to this effect, exhibiting a 
stronger suppression at low-£ The limited amount of 
available data makes it difficult to resolve these spectra 
at the present time, however, it is plausible that improved 
polarization measurements will be able to detect this ef- 
fect in the future.. 

This analysis can be easily extended to models with 
a tensor component. Since tensor fluctuations are de- 
scribed by a massless scalar field, they exactly satisfy the 
initial conditions Eq. (|43j) . Scalar and tensor modes on 
horizon and superhorizon scales at the onset of inflation 
experience the same degree of suppression, with the re- 
sult that the tensor-to-scalar ratio on these scales remains 
constant, r(k < ki—%) ~ r{kt=2). This is in marked con- 
trast to what is expected in other pre-inflation scenarios. 
Recently, Nicholson et al. showed that an initially 
fast-rolling inflaton, while yielding similar scalar spectra 
to those presented here, results in a tensor-to-scalar ratio 
that increases in magnitude on large scales. Therefore, 
while improved E-mode polarization measurements will 
be needed in order to detect pre-inflationary physics, it 
appears that knowledge of the tensor component will be 
necessary for determining what the pre-inflationary con- 
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FIG. 3: C,f E -spectra of the models presented in Figure 2. The 
pre-inflation phase has a reduced effect on the EE spectra. 



[Or 




FIG. 4: Cg -spectra of the models presented in Figure 2. The 
suppression is more pronounced in these spectra than in the 
TT- and EE-spectra, making this effect potentially observable 
with future data. 

ditions were. 



CONCLUSIONS 

We have considered the effects on the primordial power 
spectrum of a period of radiation-dominated expansion 
prior to the inflationary epoch. We impose initial con- 
ditions on the quantum fluctuations at the start of in- 
flation appropriate to the pre-inflationary era. We find 
that if the inflaton is slowly rolling during RD then it 
behaves as a massless scalar, decoupled from the pre- 
inflationary curvature perturbations. This allows for the 



construction of an adiabatic vacuum that is exact at all 
wavelengths. The largest-scale modes that exit the hori- 
zon at the start of inflation are most strongly effected 
by the phase transition, since such modes were never in 
the ultraviolet limit. This forces one to consider that 
the initially large-scale fluctuations are vacuum fluctua- 
tions of very low momentum that originate in the adi- 
abatic vacuum of the radiation-dominated era. These 
relatively low momentum modes lead to a strong sup- 
pression of the power spectrum on the largest scales. 
We emphasize that our assumption of a decoupled scalar 
evolving in a background of thermal radiation fluctua- 
tions is self-consistent, but is certainly not the only pos- 
sible set of conditions for the pre-inflationary phase. Ab- 
sent a model for the pre-inflationary universe, one could 
in principle choose initial conditions with an arbitrary 
pre-inflationary curvature perturbation. This is true for 
any inflationary model: Trodden and Vachaspati, for ex- 
ample, have shown that the universe must be homoge- 
neous on scales larger than the horizon size in the pre- 
inflationary phase in order for inflation to happen at all 

The resulting power spectrum is shown to lower 
the Cg's of the CMB temperature and polarization 
anisotropics at low £, particularly around the quadrupole. 
We find a spread in likelihoods of | — 2Aln£| < 2 in models 
for which at most 1 e-fold of inflation occurs between the 
onset of inflation and the time that scales corresponding 
to the quadrupole leave the horizon, as compared to the 
best-fit WMAP model with no tensors and no running. 
We find an improvement of only — 2Aln£ = 1 for models 
with N 60.5, so current data does not strongly favor 
such models. Since the modifications to the C^-spectra 
occur at the cosmic variance dominated low-^ multipoles, 
a more accurate measurement of the more sensitive TE 
spectrum will be necessary to decisively detect this effect. 

Hirai and Takami [lOj apply a methodology very sim- 
ilar to that presented here to investigate the effect of 
pre-inflationary physics, including that from a prior 
radiation-dominated phase. However, instead of consid- 
ering a modification to the boundary condition of the 
quantum scalar mode 6<p, they use a boundary condition 
appropriate to a hydrodynamical radiation-dominated 
mode, with c 2 s = 1/3. This results in the apparently 
unphysical conclusion that modulations to the inflation- 
ary power spectrum persist at all scales, regardless of 
the duration of inflation. However, we argue that the 
dominant contribution to the post-inflationary curvature 
perturbation will be from those generated by fluctuations 
in the inflaton field itself, for which c 2 s — 1 even in the 
pre-inflationary phase. 

Finally, we draw attention to previous work on this 
subject, namely 0] and [1, ||, who consider a pre- 
inflationary fast-rolling phase and matter dominated 
phase, respectively. While the details of the pre- 
inflationary phase differ amongst these models, the ef- 
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feet on the primordial power spectrum is strikingly sim- 
ilar. However, Nicholson et al. showed recently 
that models for which there is a transition from an ini- 
tial fast-rolling stage to the slow-roll attractor have a 
relatively large gravitational wave contribution on large 
scales, in contrast to models preceded by radiation or 
matter-dominated expansion. Our analysis differs fur- 
ther in that the suppression of power on large scales is 
manifestly the result of a choice of vacuum. The details 
of the background dynamics do not enter into the cal- 
culation, except to specify the vacuum state in the pre- 
inflationary phase. The similarity of the modifications to 
the power spectrum suggests the conclusion that the ef- 
fect of pre-inflationary physics can be generically consid- 
ered as an effect of vacuum definition, in a fashion similar 
to to that for trans-Planckian modulations [l5[ . (This is- 
sue was considered in an effective field theory context in 
Ref. |L3|.) However, in all cases the resulting modifica- 
tion to the spectrum of fluctuations in the CMB results in 
only a modest improvement of the fit to the data. Since 
the observational uncertainties at these scales are dom- 
inated by cosmic variance, evidence for these effects is 
likely to remain inconclusive in the absence of additional 
observable evidence. 

We thank Andrei Linde for helpful comments. This 
research is supported in part by National Science Foun- 
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